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Abstract 

The entanglement entropy of v = 1/2 and v = 9/2 quantum Hall states in the presence of 
short range disorder has been calculated by direct diagonalization. Spin polarized electrons are 
confined to a single Landau level and interact with long range Coulomb interaction. For v = 1/2 
the entanglement entropy is a smooth monotonic function of disorder strength. For v = 9/2 the 
entanglement entropy is non monotonic suggestive of a solid-liquid phase transition. As a model 
of the transition at v = 1/2 free fermions with disorder in 2 dimensions were studied. Numerical 
evidence suggests the entanglement entropy scales as L rather than the LlnL as in the disorder 
free case. 
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I. INTRODUCTION 



Entanglement entropy is a quantity, which, broadly speaking, quantifies the quantum 
mechanical entanglement between a system and a subsystem. It may be useful as a numer- 
ical means to detect quantum phase transitions, in particular in cases for which an order 
parameter is non local or simply unknown^. For example, for the quantum Hall state at 
filling v = 5/2, i.e. where the N — 1 Landau level is half filled, there appears to be a cusp, 
a lack of smoothness in the entanglement entropy, as a function of disorder strength^. This 
lack of smoothness, is apparent, in numerical calculations, even for small system sizes and 
is taken to be a signature of a quantum phase transition away from the Moore-Read state 
as disorder increases. 

In this paper, the behavior of the entanglement entropy as a function of disorder strength 
is investigated for other half-filled systems, in particular for filling 1/2, the N = Landau 
level being 1/2 filled, and for filling 9/2, a 1/2 filled N = 2 Landau level. Since the physics 
of the N = 0,1,2 Landau levels is quite different, quite distinct behavior as a function of 
disorder strength is anticipated. In the lowest Landau level (N = 0), the composite fermion 
picture provides a very good description of the essential physics. In the half filled case each 
electron takes two flux quantum leaving composite fermions in two dimensions moving in 
zero magnetic field 4 . Although disorder can have drastic effects in two dimensions, due to the 
limited system sizes accessible with direct diagonalization, it is anticipated there is smooth 
behavior of the entanglement entropy as a function of disorder strength. In the N = 2, third 
Landau level, the electrons form stripes and bubbles or possibly liquid crystalline states^ - -. 
In experiment, the stripes are very sensitive to disorder. (More precisely, the anisotropy 
in the diagonal resistivity is sensitive to sample quality). One therefore conjectures, in a 
numerical calculation one should see something analogous to a solid liquid transition, i.e. a 
signature of a first order phase transition. 

II. ENTANGLEMENT ENTROPY AS A FUNCTION OF DISORDER 
STRENGTH FOR v = 1/2 AND 9/2 

The numerical method used is direct diagonalization applied to square (aspect ratio one) 
clusters with periodic boundary conditions, the square torus geometry. As is standard, the 
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Landau gauge is taken for the vector potential. Spin-polarized electrons are confined to a 
single Landau level and interact with the Coulomb potential. The random potential U(r) 
is taken to be delta correlated < U(r)U(r') >= Uq5(t — r') and the disorder strength is 
given by a parameter Ur = ^3Uq/2. This approach has previously been used to study 
the entanglement entropy without a disorder potential 2 and with a disorder potential at 
fillings v = 1/3 and v = 5/2^. Due to lack of translational invariance and the necessity of 
disorder averaging the largest system size considered was 13 electrons in 26 orbitals with 
a state space of approximately 1X10 7 . Since in the third Landau level, the states are not 
necessarily liquid like, it would of interest to use clusters of different aspect ratios and pick 
the aspect ratio giving the lowest disorder averaged energy^. However, because of computer 
time limitations this was not attempted. 

The entanglement entropy is calculated by taking a subsystem of I adjacent orbitals. 
Recall in the Landau gauge, the orbitals consist of strips oriented along, say the y-axis, 
of width the order the magnetic length. The reduced density matrix is straightforward to 
calculate from the ground state wave function. It is then diagonalized, giving the eigenvalues 
Xj from which the / orbital entanglement entropy S(l) , S(l) = — J2j \ m Xj is obtained. 
This procedure is done for every realization of the random potential; the results are averaged 
to give < S(l) > where <> denotes average over the random potential. For the smallest 
systems (10 electrons in 20 orbitals) we have averaged over 200 realizations of the random 
potential; the largest systems, we have averaged over as few as 30 realizations. 

Considering first, the lowest Landau level, in figure la < 5(2) > is plotted vs. disorder 
strength for filling 1/2. 1 = 2 corresponds to a real space width for the largest system 
size (N=26) of 2^/2^/26. This is roughly one magnetic length. The system sizes 10/20 
(10 electrons in 20 orbitals), 11/22, 12/24 and 13/26 were treated. Error bars are given 
by the root mean square (rms) values of S(l), i.e. o = y ^ <s ^ l '^~_^ s ^ > with N s being the 
number of samples. One sees rather smooth behavior vs. disorder strength in comparison 
to v = 5/2, see fig 8 a of reference^. Similar behavior is obtained for larger subsystem 
sizes, in figure lb < 5(12) > is plotted for Ur for filling 1/2. For the system sizes we could 
study numerically there is not much difference in considering a subsystem of fixed I or a 
fixed aspect ratio subsystems. That is, even quantitatively, there is not much difference in 
comparing < 5(12) > for all system sizes, as we have done, or comparing < 5(10) >jv=2o, 
< 5(11) >at=22) < 5(12) >7v = 24, < 5(13) >Af=26- Once again, there is rather smooth 
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behavior of the entanglement entropy vs. disorder strength. This is at least consistent with 
what one would expect from non interacting composite fermions. In two dimensions for non 
interacting fermions, all states are exponentially localized. However, if the system size is 
smaller then the localization length, the states will appear to be extended^ 1 ^. We conjecture 
for the small system sizes and disorder strengths we study, this is the case. Experimentally^ 
it appears the composite fermi liquid is more robust to the presence of disorder then the 
quantum Hall state at filling 5/2. We also note that qualitatively, the entanglement entropy 
vs. disorder curves for v = 1/2 and v — 1/3 are similar (see figure 7 of reference^). 

Turning now to the third Landau level, figure 2a is a plot of < 5(2) > vs. disorder and 
figure 2b is a plot of < 5(12 > vs. disorder for system sizes 10/20, 11/22,12/24 and 13/26. 
For both < 5(2) > and < 5(12) > there is non monotonic behavior of the entanglement 
entropy as a function of Ur. This is reminiscent of a first order phase transition (Maxwell 
construction) which, naturally, in the third Landau level can be interpreted as a transition 
from a solid to a liquid like stated. 

III. ENTANGLEMENT ENTROPY OF FREE FERMIONS WITH DISORDER AT 
HALF FILLING 

We attempted for v — 1/2 and fixed I to perform a 1/N extrapolation (N is the system 
size). Regrettably, this was not successful due to the nonlinear behavior with 1/N. Note 
that in calculating the entanglement entropy for strips of I orbitals (length y/2/KN, width 
l^2n/N ), the subsystem, the strips, have a state space of dimension 2 l . Hence the reduced 
density matrix has a dimension 2'X2^ and the entanglement entropy is bounded from above 
by I In 3^. It is thus impossible for the entanglement entropy to scale as the length of the 
strip, as is greater then / for large N. By analogy to the area law for fermions^"—, we 
have investigated for disordered systems at 1/2 filling in the lowest Landau level, whether 
the entanglement entropy scales as Z 1//2 , I or /^lnZ. We were unable to come to any 
definite conclusion due to the small system sizes. However, it seems, not without interest 
to study a simpler problem of the entanglement entropy of free fermions in zero magnetic 
field at 1/2 filling with disorder in two dimensions. The motivation for this model is the 
composite fermion picture where electrons in the 1/2 filled Landau level are equivalent to 
non interacting composite fermions in zero magnetic field. That is, let us consider the 
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usual Anderson tight binding Hamiltonian^ii^ in two dimensions on a square lattice with 
periodic boundary conditions. There is nearest neighbor hopping t with t — 1 and a site 
diagonal matrix element that is rectangularly distributed with width W. We take the system 
to be square, with periodic boundary conditions of dimension 60 X 60 and 80 X 80 and the 
subsystem to be small squares of dimension L X L where L ranges from 1 to 30. The 
entanglement entropy was calculated with Peschel's method^ for a given realization of 
disorder and then an average over disorder is done giving < S(LxL) >. Recall, even though 
the system is half filled, since the electrons are non interacting, this only involves one electron 
calculations and hence quite large systems and subsystems can be treated. In figure 3a 
< S(LxL) > is plotted vs. LlnL for W =15. For this value of W, the localization length 
is about 1 in units of the lattice constant (see the inset of fig. 1 of ref.— ) and hence one 
anticipates both the system and subsystem to see the effect of localization. An average 
over a 1000 realizations of the disorder was taken, the error bars are barely discernible in 
the figure. Recall, although all states are localized in 2 dimensions by disorder, W=15 was 
chosen to avoid complications due to finite size effects. One notes a definite downward 
curvature in the graph 3a. This is in contrast with figure 3b in which < S(LxL) > is plotted 
vs. L again for W=15 where one sees more linear behavior. We suggest that non interacting 
states with disorder in dimension 2 obey an area law without multiplicative logarithmic 
corrections. This result is a surprise in that 2 gapless fermi systems coupled by a weak link 
have an entropy proportional to In iM- , hence multiple links should give a LlnL entropy. In 
addition in one dimension, for fermions, whether or not disorder is present, the entanglement 
entropy is proportional to In L^LZIL . 

IV. REMARKS ON THE TOPOLOGICAL ENTANGLEMENT ENTROPY FOR 
v = 5/2 WITHOUT DISORDER 

Recently, it has become possible to compute the entanglement entropy for 18 electrons in 
36 orbitals, with no disorder, in the flat torus geometry, without extreme numerical effort. 
This is due to improvements in hardware; workstations with 64 GB of RAM are not very 
costly. This technological improvement has enabled us to reconsider the following puzzling 
situation. If system sizes 24, 26,28,30 and 32 are used to extrapolate the entanglement 
entropy S(l) and then the topological entanglement entropy 7 is extracted from the S(l) 
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vs. \/l line, a value of 7 = 2.07 ± 0.15 is obtained. However, if the same procedure is 
followed including size 34, a "worse" value is found 7 = 1.93 ± 0.16^ in comparison to the 
value 7 = 2.08 for the Moore-Read state^ 1 ^. In figure 4a, the extrapolated entanglement 
entropy using system sizes 12-18 is plotted vs. l 1/2 . The error bars in figure 4a correspond 
to an estimate of an error in the linear extrapolation. The value of the y- intercept, —7 = 
—2.08 ±0.19. Thus by including system size 36 much closer agreement with the Moore-Read 
value is obtained. What is the origin of this behavior? In figure 4b, S(l) is plotted vs. l 1 ^ 2 
for system sizes 32, 34 and 36. It is hard to see system size 34 in that it is covered up 
by the system size 32 points. The growth in the entanglement entropy occurs only when 
a pair of electrons is added. This numerical effect is consistent with the Moore-Read state 
being a paired state of fermions. Returning to figure 4a, the extrapolated entanglement 
entropies obtained from even number of electrons (system size 24, 28,32,36) is plotted as 
diamonds. Fitting these points to a line, one obtains a y-intercept of —7 = —2.20 ± 0.19 
again consistent with the topological entanglement entropy of the Moore- Read state. These 
results agree with the analysis using the entanglement spectrum^— applied to the quantum 
Hall spherical geometry, that the essential physics of v = 5/2 is given by the Moore- Read 
state. 

Finally, it was possible to obtain one more system size for 7/3 filling, 13 electrons in 
39 states. In figure 5, S(l) is plotted vs. l 1/2 for various system sizes. There is little 
change in the entanglement entropy is going from 8 to 9 electrons and 12 to 13 electrons. 
Since 4 divides both 8 and 12, this behavior is consistent with the 4 fermion clusters in 
the k=4 Read-Rezayi^ state. This result, though consistent with our previous work 1 -, is 
in disagreement with recent experiments 26 that show that v — 7/3 does not appear to be 
the k=4 Read-Rezayi state. Note however, that our calculations do not include disorder; we 
expect a state with 4 fermion clusters to be very sensitive to disorder given the sensitivity 
of the Moore- Read state to disorder-^. 

V. CONCLUSIONS 

The entanglement entropy of v — 1/2 and v = 9/2 quantum Hall states in the presence 
of short range disorder has been calculated by direct diagonalization. For v = 1/2 the 
entanglement entropy is a smooth monotonic function of disorder strength. For v = 9/2 the 



6 



entanglement entropy is non monotonic suggestive of a solid-liquid phase transition. As a 
model of the transition at v — 1/2 free fermions with disorder in 2 dimensions were studied. 
Numerical evidence suggests the entanglement entropy scales as L rather than the L In L 
as in the disorder free case. It therefore may be of some interest to use the entanglement 
entropy as a numerical tool to study interacting electrons in zero magnetic field but in the 
presence of disorder, particularly in two dimensions. 

We thank Mahito Kohmoto for a helpful conversation. This work was supported in part 
by NSF Grant no. 0705048 (B.F. and D. L.) and the Department of Energy DE-FG02- 
08ER64623— Hofstra University Center for Condensed Matter (G. L.). 

Figure Captions 

figure 1. < S(l) > vs. Ur for filling 1/2. In figure la < 5(2) > is plotted, while figure lb is 
a plot of < 5(12) >. 

figure 2. < S(l) > vs. Ur for filling 9/2. In figure 2a < 5(2) > is plotted, while figure 2b is 
a plot of < 5(12) >. 

figure 3. < S(LxL) > for an Anderson model on a two dimensional square lattice at 1/2 
filling. W = 15 and the system sizes are 60 x 60 and 80 x 80 with periodic boundary 
conditions. In figure 3a S(LxL) vs. L In L is plotted. In figure 3b S(LxL) vs. L is plotted. 
An average over a 1000 realizations of the disorder was taken, the error bars are barely 
discernible in the figures. 

figure 4. Entanglement entropy vs. I 1 / 2 for v = 5/2, no disorder. In figure 4a the extrapo- 
lated entanglement entropy is plotted. In figure 4b the entanglement entropy for system sizes 
16/32, 17/34 and 18/36 is plotted. The error bars in figure 4a correspond to an estimate of 
an error in the linear extrapolation. 

figure 5. 5(/) vs. I 1 / 2 for various system sizes v = 7/3. 
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figure 4a 
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figure 4b 
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figure 5 
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